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We study in detail the structure of phase space in the neighborhood of stable periodic orbits in a 
rotating 3D potential of galactic type. We have used the color and rotation method to investigate 
the properties of the invariant tori in the 4D spaces of section. We compare our results with 
those of previous works and we describe the morphology of the rotational, as well as of the tube 
tori in the 4D space. We find sticky chaotic orbits in the immediate neighborhood of sets of 
invariant tori surrounding 3D stable periodic orbits. Particularly useful for galactic dynamics 
is the behavior of chaotic orbits trapped for long time between 4D invariant tori. We find that 
they support during this time the same structure as the quasi-periodic orbits around the stable 
periodic orbits, contributing however to a local increase of the dispersion of velocities. Finally 
we find that the tube tori do not appear in the 3D projections of the spaces of section in the 
axisymmetric Hamiltonian we examined. 

Keywords: Chaos and Dynamical Systems, Galactic Dynamics, 4D surfaces of section, Invariant 
tori 
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1. Introduction 

The method of surfaces of section for the study of 
dynamical systems dates back to Poincare (1892) 
and has many applications to Dynamical Astron- 
omy (for a review see e.g. Contopoulos 2002). A 
basic problem in Hamiltonian Systems of three de- 
grees of freedom is the visualization of the sur- 
faces of section. Let us assume the phase space of an 
autonomous Hamiltonian system, that has 6 dimen- 
sions, e.g. in Cartesian coordinates, (x,y, z,x,y, z). 
For a given value of the Hamiltonian a trajectory 
lies on a 5D manifold. In this manifold a surface of 
section is 4D. This does not allow us to visualize 
it directly. Several methods have been applied to 
overcome this problem in the past and we summa- 
rize them below. 

The structure of the 4D space phase space was 
examined for the first time in the pioneer work of 
Froeschle (1970, 1972). In that work he used stereo- 
scopic views and the method of slices in order to 
understand the structure of tori, that appeared at 
the neighborhood of stable periodic orbits. Similar 
methods have been applied by Martinet and Magne- 
nat (1981), Contopoulos et al. (1982), and Magne- 
nat (1982) for studying the 3D projections of invari- 
ant tori in the 4D surface of section or in the phase 
space of a 4D symplectic map. The 2D projections 
of such invariant tori have been examined on var- 
ious 2D planes in detail (e.g. Skokos, Contopoulos 
and Polymilis 1997, 1999). 

In the present paper we use the method of 
color and rotation proposed by Patsis and Zachi- 
las (1994). In this method we first consider 3D pro- 
jections and rotate the 3D figures on a computer 
screen to observe the figure from all its sides. Then 
we use colors to indicate the 4th dimension. For this 
purpose we make use of the "Mathematica" package 
(Wolfram 1999) and of its subroutines. Each point is 
colored according to the value of its 4th coordinate 
in the following way: Firstly we define the surface 
of section that we will use, e.g. y = with y > 0. 
Secondly we select a 3D subspace of the surface of 
section, e.g. (x,x,z) and we determine the mini- 
mum and maximum values of the 4th coordinate z. 
Finally we normalize the resulting interval [min(z), 
max(z)] into [0,1] from which interval the color func- 
tions of "Mathematica" take values. In our figures 
we give always the color function of "Mathemat- 



ica" that we used in a color-coded bar. The point of 
view of the observer of the figures is given in spher- 
ical coordinates. This specifies the direction along 
which the figure is projected. The unit for distance 
d of the consequents of the surface of section from 
the observer is given by "Mathematica" in a spe- 
cial scaled coordinate system in which the longest 
side of the bounding box, which surrounds the fig- 
ure, has length 1. For all figures we use d = 1. The 
method associates the smooth distribution or the 
mixing of colors, with specific types of dynamical 
behavior in the 4th dimension (Patsis and Zachilas 
1994). 

In order to study the structure of phase space at 
the neighborhood of a periodic orbit (p.o.), we first 
locate it by means of an iterative method and cal- 
culate its stability. The calculation of the linear sta- 
bility of a periodic orbit is based on the method of 
Broucke (1969) and Hadjidemetriou (1975). We first 
consider small deviations from its initial conditions 
and then integrate the orbit again to the next up- 
ward intersection. In this way a 4D map (Poincare 
map) is established, which is unique (Abraham and 
Marsden 1978 p. 521) and relates the initial with 
the final point. The relation of the final deviations 
of this neighboring orbit from the periodic one, with 
the initially introduced deviations can be written in 
vector form as £ = M £o- Here £ is the final devia- 
tion, £o is the initial deviation M is a 4 x 4 matrix, 
called the monodromy matrix. This matrix satis- 
fies the symplectic identity and the map is called 
symplectic (Arnold and Givental 2000). It can be 
shown that the characteristic equation can be writ- 
ten in the form A 4 + aA 3 + (3\ 2 + aA + 1 = 0. Its 
solutions A,, i = 1,2,3,4, due to the symplectic 
identity of the monodromy matrix, that obey the 
relations A1A2 = 1 and A3A4 = 1 can be written as 



1 -fc^y^i 

Al 'AT" 2 

. 1 -b 2 ± ^M^l 

A3 'A^" 2 



where 



a± y/A 



'1,2 



(1) 



(2) 



and 



1 throughout the paper we will refer to "n-dimensional" as 
"nD"; i.e. 3D, 4D etc. 



A = a 2 - 4(/3 - 2) 



(3) 
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The quantities b\ and 6 2 are called the stability 
indices. Following the notation of Contopoulos and 
Magnenat (1985) if A > 0, \b x \ < 2 and |6 2 | < 2, all 
four eigenvalues are complex on the unit circle and 
the periodic orbit is called "stable" (S). If A > 
and I fox I > 2, |6 2 | < 2 or |&i| < 2, |6 2 | > 2, the 
periodic orbit is called "simple unstable" (U). In 
this case two eigenvalues are on the real axis and 
two are complex on the unit circle. If A > and 
I b\ I > 2 and |6 2 | > 2, the periodic orbit is called 
"double unstable" (DU) and the four eigenvalues 
are on the real axis. Finally if A < the peri- 
odic orbit is called "complex unstable" (A). In this 
case the four eigenvalues are complex numbers and 
they are off the unit circle. For the generalization of 
this kind of instability in Hamiltonian systems of N 
degrees of freedom the reader may refer to Skokos 
(2001). When two eigenvalues collide at (1,0) of the 
unit circle the parent family becomes simple unsta- 
ble (U) and a new family of p.o. is born. In this 
paper we examine the evolution of the phase space 
at the transition from stability to simple instability 
(S — > U). The parent family in our example has an 
orbital plane, i.e. it is 2D, and by becoming sim- 
ple unstable as the energy increases, it generates by 
bifurcation a stable 3D family of p.o.. 

According to the KAM theorem (Kolmogorov 
1954, Moser 1962, Arnold 1963) in an almost In- 
tegrable Hamiltonian system of N degrees of free- 
dom there are orbits that lie on N-dimensional 
tori. The KAM theorem has been extended to 2n- 
dimensional almost integrable symplectic maps by 
Wiggins (2003, p. 225) and by Kuksin and Poschel 
(1994). This means that in a 2n-dimensional almost 
integrable symplectic map there are orbits that lie 
on n-dimensional tori. In our 3D Hamiltonian sys- 
tem a 4D Poincare map is defined on the surface of 
section. In specific cases this 4D Poincare map is al- 
most integrable. Then, according to Wiggins (2003, 
p. 225), we have orbits that lie on 2D invariant tori 
in the 4D space of section. The structure of these 
tori, is the subject of our present, rather descriptive, 
paper. 

In order to facilitate the discussion of the figures 
in the paper, we give some useful definitions related 
to a torus. In the 3D space, a torus is a surface that 
is generated when we rotate a circle on the x-z (or 
y-z) plane around the axis z (Fig. [T]). The definition 
of the angles u and v are also noted in Fig. [TJ The 
internal surface of the torus is defined as the set of 
points of the torus, where we have 90° < u < 270°. 
The external surface of torus is defined as the set of 



points of the torus, where we have 0° < u < 90° and 
270° < u < 360°. We will call "tori" the objects we 
study in this paper, despite the fact that they seem 
to be generated by an ellipse instead of a circle i.e. 
they are elliptic tori. 



z 




y 



Fig. 1. The definition of the v and u angles on a torus. 
Arrows indicate the direction along which the angles u,v in- 
crease from zero towards larger values. 



The purpose of this paper is to study the struc- 
ture of the phase space in the neighborhood of sta- 
ble periodic orbits in 3D galactic potentials. We 
want to understand how the phase space is struc- 
tured and how it evolves as the stability of the xl, 
2D central family of our system on the equatorial 
plane, changes from stable to simple unstable. We 
want also to understand how the phase space struc- 
ture evolves as the main parameters of our system 
vary. In section 2 we describe our Hamiltonian sys- 
tem, while in section 3 we present the morpholog- 
ical evolution of the families of periodic orbits in- 
volved in our study. Section 4 discusses the struc- 
ture of the phase space when we perturb the initial 
conditions of our stable periodic orbits in the x- 
direction. We study successively the dynamical be- 
havior close to the periodic orbit before (Sect. 4.1) 
and after (Sect. 4.2 ) the transition of the central 
family from stability to instability. We discuss as 
well the "rotation numbers" we define on the tori 
we have found (Sect. 4.3). Then, we increase the en- 
ergy and we describe the changes we observe in the 
spaces of section in Sect. 4,4. The perturbations in 
the z-direction are presented in Sect. 5, which has a 
similar structure as Sect. 4. In Sect. 6 we vary the 
perturbation of our system. Finally in section 7 we 
discuss our results and we enumerate our conclu- 
sions. 
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2. The Hamiltonian System 

The system we use for our applications rotates 
around its z-axis with angular velocity f2b . The 
Hamiltonian of the system in Cartesian coordinates 
is : 

H(x,y,z,x,y,z) = 
^(x 2 + y 2 + z 2 ) + $(x,y,z) 

-\^ 2 + y 2 ) (4) 

where $(x, y, z) is the potential we used in our 
applications, i.e.: 

®(x,y,z) = 
GMi 

~(* 2 + 4 + [a 1 + (4 + &2)l/2 ] 2 ) l/2 " 

GM 2 

(x 2 + $ + [a2 + (4 + qy/^y/i 

la % 

(5) 

The potential in it's axisymmetric form (q a = 
l,qt) = 1) can de considered as a representa- 
tion of the potential for the Milky Way approx- 
imated by two Miyamoto disks with masses M± 
and M 2 respectively (Miyamoto and Nagai 1975). 
In our units, distance R=l corresponds to 1 kpc. 
For the Jacobi constant (hereafter called the "en- 
ergy") Ej=l corresponds to 43950 (km/ sec) 2 . We 
have used the following values for the parameters: 
01 = kpc, hi = 0.495 kpc, M\ = 2.05 x 
10 10 M Q , a 2 = 7.258 kpc, b 2 = 0.520 kpc, M 2 = 
25.47 x 10 10 M Q , q a = 1.2, q b = 0.9 and Q b = 
60 km s~ x kpc~ 1 . The parameters q a ,qb determine 
the geometry of the disks, while a, b are scaling fac- 
tors (Binney & Tremaine 2008, p. 73- 74). The chosen 
fi& value puts corotation at R=4/cpc and has been 
used by Englmaier and Gerhard (1999). 

In Fig. [2] we give the (Ej,x) Zero Velocity Curve 
"ZVC" (see e.g. Contopoulos 2002 p.391). It refers 
to orbits on the equatorial plane z=0 and sepa- 
rates the regions where motion is allowed from those 
where it is forbidden. 

3. The orbital evolution along a 
S — > U transition 

A method to follow the stability of a family of peri- 
odic orbits in a system is by means of the " stability 
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Ej 



Fig. 2. The (Ej,x) Zero Velocity Curve (ZVC), for orbits 
on the z =0 plane, in our Hamiltonian system (4), with $ 
as in (5) and parameter values given in the text. Motion is 
forbidden in the grey areas. L4 and L5 are the Lagrangian 
points on the minor axis of the galaxy. 
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Fig. 3. Stability diagram for —5.5 < Ej < —4.5, that shows 
the stability of the family xl (black lines for its bl,b2 indices) 
and its bifurcating families of p.o. xlvl (red) and xlv2 (blue). 
The two indices of the family xlvl join at Ej =— 4.62 and the 
family becomes complex unstable. 

diagram" (Contopoulos & Barbanis 1985; Pfenniger 
1985a). The stability diagram gives the evolution of 
the stability of a family of periodic orbits in a sys- 
tem as one parameter varies, by means of the evo- 
lution of the stability indices bl, b2. In our case the 
parameter that varies is the energy Ej . Fig. [3] gives 
the evolution of the stability of the central family of 
periodic orbits in our system, xl (Contopoulos and 
Papayannopoulos 1980), and its bifurcations at the 
interval —5.5 <Ej< —4.5. We observe that xl (black 
lines for its indices bl and b2) is initially stable and 
at Ej=— 5.1644 it becomes simple unstable. There 
we have a S — > U transition and a new family, xlvl 
(red lines), is bifurcated and is stable. We call the 
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Fig. 4. Periodic orbits of xl at a) Ej= -5.624726 (before 
"A"), b) Ej= -5.010526 (between "A" and "B") and c) Ej= 
-4.732626 (beyond "B"). 




0.5 0.5 1 1.5 



Fig. 5. A typical orbit of xlvl at Ej = -5.125377 (between 
"A" and "B"). The (x,z) and (y, z) projections are given 
in enlarged scale, in order to better view the corresponding 
morphology. 



transition point "A". The family xl becomes sta- 
ble again at Ej = —4.98, point "B", where we have 
an U — > S transition for xl. At "B", is introduced 
the family xlv2 (blue line), which is initially simple 
unstable. 

In Fig. |4] we present the morphological evolu- 
tion of xl as it appears before "A", between "A" 
and "B" , and beyond "B" . Point "A" is associated 
with the vertical 2/1 resonance (for a definition see 
e.g. Contopoulos 2002 p. 379). The importance of 
the 3D bifurcating family introduced at the vertical 
2/1 resonance for the structure of galactic disks has 
been underlined in several studies (Pfenniger 1985b; 
Patsis and Grosb0l 1996; Skokos et al 2002a,b; Pat- 
sis et al 2002). Here we adopt the nomenclature 
of Skokos et al 2002a, b and we call the bifurcating 
families at Ej = —5.1644 and Ej = —4.98 xlvl and 
xlv2 respectively. In Fig. [5] we give the morphology 
of xlvl in a typical orbit at energies between "A" 
and "B". 

4. Perturbations parallel to the 
equatorial plane 

4.1. Spaces of section before the 
S — >■ U transition 

First we examine the surfaces of section for our sys- 
tem at energies before "A". Figs. and [6b de- 
scribe the surface of section for Ej=— 5.207 (before 
"A"). Despite the fact that the system we investi- 



gate is of galactic type, the main goal of the present 
study is to understand the structure of the phase 
space in various cases in a 3D autonomous Hamil- 
tonian system. In that sense we integrate our orbits 
for times necessary to obtain a clear view of the 
dynamical phenomenon we study, regardless of the 
physical meaning of the integration time interval, 
which can be more than a Hubble time. In Fig. [6^ 
we observe the invariant curves around periodic or- 
bits of two 2D families located on the equatorial 
plane z=0, in the (x, x) space with initial conditions 
(s + Ax ,x ,z ,z ) = (0.18312784 + Ax , 0,0,0) 
with Axo = 0.1, . . . , 0.8 and (xq + Axo, xq, zq, zq) = 
(-0.59595941 + Ax, 0,0,0) with Ax = 0.1,... ,0.4. 
Every invariant curve consists of 10 3 consequents. 
The diagram describes a typical situation in ro- 
tating galactic potentials. The stable periodic or- 
bit for x > belongs to xl (initial conditions 
(x , x , z , zq) = (0.18312784, 0, 0, 0)) and the stable 
periodic for x < to the retrograde family x4 (see 
Contopoulos 2002, p. 391), which has initial con- 
ditions (x ,x ,z ,io) = (-0.59595941,0,0,0). The 
extent of the invariant curves is limited by the ZVC. 
Varying the initial condition xo above the upper and 
below the lower limit of the ZVC at this Ej, always 
considering the plane of section y = (cf . Fig. [2]) , 
we find the expected bell- type curves (10 3 points 
in the surface of section), which are related with 
escape orbits, as found by Contopoulos & Patsis 
(2006) (Fig. [6b). 

4.2. Spaces of section after the 
S — >■ U transition 

Proceeding beyond "A" , between "A" and "B" , e.g. 
at Ej=— 5.1574, we encounter, always for the surface 
of section y = 0, y > 0, two simple periodic or- 
bits with positive xq. They are xl (U), with initial 
conditions (xq,xq, z,zq) =(0.18958522,0,0,0) and 
xlvl (S), with initial conditions (xo, xo, Zq, Zq) = 
(0.18939859,0,0.030236585,0). We investigate the 
phase space structure close to these two periodic 
orbits, firstly by perturbing only the x initial con- 
ditions of xl and xlvl by Ax = 0.1, 0.2 . . . 0.7 suc- 
cessively. The 2D simple unstable periodic orbit xl 
lies on the equatorial (x, y) plane of our galactic 
model. By perturbing the initial conditions only in 
x and keeping the rest equal to 0, we encounter 
non-periodic orbits that remain on the (x, y) plane. 
Nevertheless, the surface of section of our 3D au- 
tonomous Hamiltonian system is 4D, (xo, xo, Zq, zq), 
and we can consider the (x,x, z) projection. For 
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Fig. 6. The (x,x) surface of section for Ej=— 5.207. (a) The 
invariant curves are around the p.o. xl (positive x) and x4 
(negative x). (b) The surface of section considering initial 
conditions with x above the upper and below the lower limit 
of the ZVC at this Ej. We observe the appearance of "bell- 
type" curves (Contopoulos & Patsis 2006). 

the simple unstable (U) periodic orbit xl, we have 
z = always, thus the [x, x, z) is identical with the 
(x, x) projection. This projection can be observed in 
Fig. where we have seven invariant curves, sur- 
rounding the fixed point with the initial conditions 
of xl. Each invariant curve has about 10 3 conse- 
quents. We name these curves Si a , S2 a , ■ ■ ■ , Sj a for 
the (xo + Axo) perturbation of the xo initial condi- 
tion with Axq = 0.1, 0.2, . . . 0.7 respectively. 

In a way this is a counterintuitive result, since 
at this Ej xl is characterized as simple unstabl^, 
having two of its eigenvalues on the real axis. We 
have also calculated the Lyapunov Characteristic 



2 Broucke (1969) characterizes this type of instability as 
"semi-instability", subdividing it to even- and odd- semi- 
instability, depending on whether the eigenvalues on the real 
axis are positive or negative respectively. In our case the 
eigenvalues are positive (even- semi-instability). 



Number (LCN) and we found it positive, equal 
to 0.025958. The Lyapunov Characteristic Number 
(LCN@ of the periodic orbits is defined as the max- 
imum <7j, where: a\ = ^ln(|Aj|), r is the period 
of the periodic orbit and Aj the eigenvalues of the 
monodromy matrix of the Poincare map (e.g. Licht- 
enberg and Lieberman 1992 p. 302, Skokos 2010). 
Figure [3 underlines the fact, that the dynamical 
behavior close to a simple unstable orbit in a 3D 
Hamiltonian system can be similar to that at the 
neighborhood of a stable periodic orbit if the per- 
turbation is restricted in one direction (in our case 
it is the radial one). 

Applying the same seven perturbations to the 
xq initial conditions of xlvl we encounter in the 
(x,x,z) space seven tori. The projection of the fig- 
ure is given in Fig. [8l For this projection we use the 
point of view which is determined in spherical coor- 
dinates by the angled (6,<j>) = (17°, 20°) (see also 
introduction). We observe seven tori drawn with 
red color surrounding the xlvl periodic orbit and 
we name them Si, S 2 , S3, S4, S5, Sq and Si starting 
with the closest to the periodic one. We also plot 
with green color the Si a , S 2a , S 3a , S^a, S 5a , 5 6a and 
Sia invariant curves. Each torus consists of 10 4 con- 
sequents. Being just after the bifurcation point, the 
x initial conditions of xl and xlvl are almost iden- 
tical. Thus, both the Si tori and the Si a invariant 
curves surround the "common" initial x value. On 
the [x, x) projection the Si tori and the Si a invari- 
ant curves practically overlap. 

The Si tori are 3D projections of the invariant 
tori of our 4D Poincare map. Their structure, except 
that of 54, resembles the morphology of the objects 
defined as "rotational tori" by Vrahatis et al. (1996, 
1997). These authors found the same kind of tori in 
a 4D symplectic map related with the problem of 
beam stability in circular particle accelerators. Even 
S*2, which seems to have its own internal structure 
(Fig. [8]), becomes a typical rotational torus if we 
integrate for time giving 4 x 10 5 consequents (Fig. 
[9]). An exception from this morphology is the torus 
54, which has a thin, complicated, ribbon-like struc- 
ture. Fig. [TUl presents S4 in the 3D (x,x,z) space. 
The figure helps us understand that 6*4 intersects 
itself at five places, namely A,B,C,D and E. In prac- 



When we use the term Lyapunov Characteristic Number we 
mean the maximal Lyapunov Characteristic Number. 
4 The (0°, 0°) projection brings the x-axis and y-axis horizon- 
tally and vertically on the plane of the paper respectively and 
the z-axis perpendicular to them. 




Fig. 8. The (x,x,z) 3D projection of the (x,x,z,z) 4D surface of section at Ej=— 5.1574. Our point of view is (6, <f>) = 
(17°, 20°). The Si tori around xlvl (see text) are drawn with red color, while the Si a invariant curves around the xl periodic 
are given with green lines. 
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Fig. 9. The (x, x, z) projection of torus S2 in the (a;, x, z, z) 
4D surface of section for Ej =— 5.1574. Our point of view in 
spherical coordinates is now (6,(j>) = (40°, 15°). This is a 
typical morphology of a 3D projection of a rotational torus. 
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Fig. 10. 3D projection of torus 54 from the point of view 
(6, 4>) = (40°, 264°). A,B,C,D and E are the five intersections 
of S4 by itself. 



tice this can be realized only by rotating the figure 
on the screen of our computer to understand its 
detailed morphology. The structure of our S4 tori 
resembles the morphology of the objects defined as 
"tube tori" again by Vrahatis et al. (1997). 

Let us have now a closer look at the 4D struc- 
ture of the rotational Si tori by applying the method 
of color and rotation. The tori occupy a subspace 
of the 4D space of section bounded by [xi,ac2] x 
[±i,x 2 ] x [zi,z 2 ] x [zi,z 2 ] = [-0.154927,0.890256] x 
[-1.79387,1.77590] x [-0.0302, 0.0302] x 
[—0.096,0.096]. As an example of the 4D structure 
of a rotational torus we depict in Fig. [TT] the first 
4 x 10 5 consequents of the 52 torus. This torus has 



a small thickness, as it is generated by rotation 
around the x-axis of a thin ellipse. We have cho- 
sen to plot the consequents in the (x, z, z) projec- 
tion and color them according to their values in the 
x-dimension. With the help of a graphic software 
we have rotated S 2 in our computer screen in order 
to view it from all different perspectives and better 
understand its internal structure. We observed that 
there is no mixing of colors on the surface of S 2 . 
On the contrary, we find that the color variation 
on it follows some rules. This property of the "rota- 
tional tori" was already known by Patsis & Zachilas 
(1994). However, the details, presented here for the 
first time, indicate a generic behavior that charac- 
terizes this class of objects. 

Moving along the u-direction of the torus for 
constant u (see Fig. [T]) the consequents keep their 
color by changing from the external to the internal 
side of the torus and vice versa when they reach 
certain four lines along the directions labeled with 
A, B, C and D in Fig. QU 

The details of the color variation close to these 
lines as we move along the v- direction for constant 
u, is shown in Fig. [12] This is an enlargement of the 
right part of Fig. [11] where e.g. red colored conse- 
quents on the external side of the torus change side 
sliding to the internal one, along the line indicated 
with "A" (upper left side in Fig. [12]). If we continue 
moving clockwise on the f-direction the red conse- 
quents, now on the internal side of the torus, con- 
tinue until we reach area "D". There, they change 
again side and turn to the external side. This al- 
lows us very clearly to observe the way that the 
transition from one side to the other happens. At 
the lower side of the figure we can observe how the 
red points slide behind the blue, forming a net. The 
"red net" is filled with red if we continue the cal- 
culation for longer time. Similarly the "blue net" 
on the inner side of the torus will be filled by blue 
after a long time and so do the red consequents at 
the lower left part of Fig. [12] will be covered by the 
blue. In our example we started moving along the 
-u-direction from a red point. In general, the suc- 
cession of colors in the t;-direction with constant u 
angles can be understood by looking at the color- 
coded bar given at the right of the figures with the 
colored tori, e.g. in Figilll At the change of side 
of colors, red is combined with violet-blue (above 
or below it), orange with blue, yellow with light- 
blue etc. and finally at u = 90° and u = 270° the 
consequents have the shade of green (middle of the 
color-code bar). 
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Fig. 11. The torus 5*2 in the (x,x,z,i) 4D surface of section for Ej=— 5.1574. The location of the consequents is given in 
the (x, z, i) projection and are colored according to their value in the x coordinate. Our point of view in spherical coordinates 
is (6,<f>) = (180°,9°). 



A 




Fig. 12. The right part of the torus in Fig. [TT]in enlargement. We observe how consequents of a certain color change from 
one side of the torus to the other. 



Besides the color variation in the ■u-direction, 
there is also color variation along u (constant v). If 
we examine e.g. a "red-dominated" area on the S2 
torus in Fig.[TT]we observe, that we get a red (x = 0) 
in the middle of the torus (u = 0). At first for 
u : 0° — > 90° we observe a color variation from red 
to yellow and then to green, for values A = — )■ 1 
(Fig. HU). Secondly for u : 90° -> 180°, at the inter- 



nal side of the torus, the colors change from green 
to light blue, then to blue and finally to violet (now 
for values A = 1 -)• 0). Then, for u : 180° -> 270° 
the colors change from violet to green (for values 
A = -1) and finally for u : 270° -)■ 360°, at the 
external side of the torus, we can observe a smooth 
color succession from green to yellow and then to 
red for values x = — 1 — >• 0. In conclusion we observe 
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a cyclic variation of colors (values of x) along the 
u-direction. This means that the color variation is 
also smooth along described above. The color 
variation is similar for all rotational Si tori. 

In other cases of our system we find elliptic tori 
with considerable smaller ellipticity. We give an ex- 
ample not belonging to the S —> U transition we use 
in our presentation. By perturbing the x initial con- 
dition of a stable periodic orbit with initial condi- 
tions (x ,i ,z,z ) =(0.6765982,0,0.254816,0) for 
Ej=— 4.3 (family xlv3, Skokos et al. 2002a, b) with 
Ax = 10 -5 we observe a torus which we name S a 
(Fig. [T3|) . In Figs. [TBI and [T41 torus S a is depicted in 
the 3D subspace (x,z, z) of the 4D surface of sec- 
tion and is colored according to the values of the 
4th dimension x. In Fig. [13] we observe that at the 
upper region of the external surface of the S a dom- 
inate blue shades. Along the u-direction we observe 
a smooth color variation from blue to light blue, to 
green and to yellow. At a region indicated with an 
arrow labeled "A", we observe that yellow meets 
the blue. At this region we have an intersection of 
the external surface with the internal surface of the 
torus S a . The details of the 4D structure of the torus 
close to the intersection is described in Fig. [T4l The 
color continues from a yellow shade on the exter- 
nal surface (indicated with "Ext") of the torus S a 
to orange and then to red on the internal surface 
of the torus (indicated with "Int" in Fig. [j~4"]) . We 
observe a smooth variation of colors from red to the 
orange on the internal surface of the torus S a (see 
also Fig. [T3|) until a region where orange meets the 
blue. At this region, indicated with arrow "B", we 
have in the 3D (x,z,z) projection again an inter- 
section of the external surface of the S a with the 
internal surface of the torus S a along a line. 

We want to apply now the color and rotation 
method to the other kind of tori we found around 
xlvl in Fig. [HJ As we already remarked, S4 has a dif- 
ferent morphology than the rest of the Si's, i.e. it is 
a tube torus. In the 3D projection (x, x, z) (Fig. [TO]) 
we have realized that £4, intersects itself in five re- 
gions (A, B, C, D and E). In Fig. [15] the £4 torus is 
colored according to the z values. In this figure we 
observe a smooth color variation at the areas be- 
tween two successive intersections. For example we 
examine the area between the intersections A and B 
on the internal surface of the torus along the arrow 
"1". If we move from A as indicated with the ar- 
row"!" (Fig. [To"]) , for constant u, and starting with 
blue colored points, we can see a smooth color vari- 
ation from blue to light blue, to green, to yellow, to 



orange at the intersection B. Now if we move from 
the intersection A along the direction indicated with 
the second arrow, "2" we have drawn at the "cross- 
road" A in Fig. [T5l we observe that the succession 
of colors starting from red, at constant u, is red — > 
orange — > yellow — > green — > light blue — > blue at 
the intersection B. In all cases we observe at the in- 
tersections of 54 in the space (x, x, z) two different 
tubes, coming from different directions and having 
different colors. For example at the intersections A 
red meets blue and at the intersection B blue meets 
orange. Between two intersections of S4 we observe 
a smooth color variation and in all cases at the in- 
tersections meet different colors. This means that 
the points have different values at the 4th dimen- 
sion and the intersections are not intersections in 
the 4D space but only in the 3D projections. 

4.3. Rotation Numbers 

A point that we want to investigate, is whether some 
properties of the rotational and tube tori are re- 
flected in quantities, that could be defined in corre- 
spondence with the rotation numbers on the usual 
invariant curves around stable periodic orbits of 2D 
systems. Thus we have first to define a rotation 
number in the 2D and 3D projections of a 4D torus. 
If our 4D space of section is (x, z, x, i), with the p.o. 
at (xo,zq,xq,zo) we can consider e.g. its 2D projec- 
tion in (x,x) and its 3D projection in (x,x,z). Our 
motivation for this is the observation, that succes- 
sive consequents form an invariant torus by filling 
its surface in different ways. Following the formation 
of an invariant torus on a screen as the number of 
consequents increases, one gets the impression that 
wires are wrapped around the surface of the torus, 
following different patterns for different tori. This 
seems to be a straightforward counterpart of the 
different angles at which an invariant curve is filled 
by successive consequents in 2D systems, which is 
described by the "rotation number" of an invariant 
curve. We want to attribute a similar number to 
each invariant torus. Let us assume that the projec- 
tion of the torus on the (x, x) plane gives an "invari- 
ant curve" , in the same way that the Sj tori appear 
projected in Fig. [7] In such a diagram the rota- 
tion number (rot) is defined as for a usual invariant 
curve around a stable periodic orbit at (xo,x'o) on 
a 2D surface of section. In fact ll rot" is the average 
rotation angle along the invariant curve. This quan- 
tity is different for different invariant curves and the 
variation of the rotation number as a function of the 
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Fig. 13. The torus S a around a stable p.o. xlv3 in the (x,x, z, z) surface of section for Ej=— 4.3. The consequents are colored 
according to their value in the x coordinate. Our point of view in spherical coordinates is (0, cj>) — (30°, 30°). We observe again 
a smooth color variation on the surface of the torus. The color changes side along the lines indicated with arrows. 
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Fig. 14. The region of the intersection of the external surface with the internal surface of the torus S a in the (x,x,z,z) 4D 
surface of section for Ej=— 4.3. The consequents are colored according to their value in the x coordinate. Our point of view in 
spherical coordinates is (9, </>) = (30°, 30°). Arrows indicate the external and internal sides of the torus. 



distance in the direction of a coordinate, e.g. the x- 
coordinate, defines a "rotation" curve (Contopoulos 
2002, pg. 139). 

Accordingly, for the calculation of rot for the 
torus in the (x, z, x) projection, 



(1) we define an initial rotation angle fo for the 
initial point of an orbit on the invariant torus 
with coordinates (xi,&i, z±) as the angle formed 
by the vector joining this point with the central 
periodic orbit (xq,xq,zo) with the plane (x,x). 
More precisely: 



fo = arctan 



z\ - z 



y/(xi - x ) 2 + (ii - x ) 2 



(6) 



(2) We compute the rotation angle fj 

between two successive consequents. The i-th 
consequent with coordinates and the 

i+l-th with coordinates (xi+i, ii+i, Zi+i) on an 
invariant torus. We define the vectors X and 

Y as X = (x i+ i - xq, x i+ i - x , z i+1 - z ) and 

Y = (xj — xq, Xi — x , Z{ — zq). respectively We 
find the rotation angle from the inner product 
of X and Y by means of the formula: 
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Fig. 15. The (x,x,z,z) 4D surface of section of torus S4 for Ej=— 5.1574. The consequents are colored according to their 
value in the i coordinate. Our view angles are (9, <j>) — (30°, 45°). 
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Fig. 16. Rotation curves at Ej = —5.1574 for the invariant curves around the periodic orbit of the xl in the 2D subspace 
(x,x) (green line), around the periodic orbit of the xlvl in the 2D subspace (x,x) (blue line), and around the periodic orbit of 
the xlvl in the 3D subspace (x, x, z) (red line). We observe that the three curves have a very similar variation, and practically 
overlap. 



h = arccos — — (7) 

(3) we finally compute the average of all fi's to get 
rot for the torus. 

The rotation curves along the x-direction, 
for the invariant curves around the simple 
unstable periodic orbit of the xl family at 
Ej=— 5.1574, with initial conditions (xq,zo,xq,zq) 
= (0.189585220,0,0,0), is given in Fig. [H (green 
curve). Notice that xl is stable as regards devia- 
tions on the plane of symmetry, thus we can define 



a rotation number for orbits close to xl, starting 
on this plane. On the x-axis of Fig. [16] we give the 
distance from the periodic orbit in the x-direction. 
Similar curves are calculated for the tori around the 
stable periodic orbit of the xlvl family at the same 
energy (Ej = —5.1574), which has initial conditions 
(x ,Zo,x ,Zo)= (0.18939859,0.030236585,0,0). Two 
rotation curves for xlvl are given also in Figs. [16] 
using the (x, x) (blue line) and (x, x, z) (red line) 
projections respectively. For the calculation of the 
rotation numbers in these rotation curves we fol- 
lowed the definitions mentioned above. We observe 
that in all cases we have a similar variation of the 
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rotation numbers. The behavior of rotation curves 
for the invariant tori around the stable periodic or- 
bit of the xlvl family has similar behavior with the 
rotation curve for the invariant curves around the 
simple unstable periodic orbit of the xl family. We 
realize that the rotation numbers of the tube tori, 
follow the rotation curve and occupy the expected 
position in this diagram without any kind of excep- 
tional behavior. 
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Fig. 17. The invariant curves T la , T 2a , T 4a , T 5a , T 6a and T 7a 
in the (x,x) surface of section for Ej=— 5.131377. 




Fig. 18. The 3D (x, z, x) projection of the 4D surface of 
section for Ej =— 5.131377. Our point of view in spherical co- 
ordinates is (8, 4>) = (18°, 23°). The invariant tori around the 
initial condition of the stable p.o. xlvl are colored red and 
the invariant curves around the simple unstable p.o. of xl on 
the z = plane are colored green. The initial conditions of xl 
and xlvl are given with green and red crosses respectively. 



4.4. Energy variation 

Coming back to the evolution of the stability of 
the families xl and xlvl, as described in Figj3l we 
study the changes introduced in the structure of 
the phase space as the energy increases. At the bi- 
furcating point the initial conditions of the fami- 
lies xl (parent family) and xlvl (bifurcating fam- 
ily) are identical. The xlvl family comes in the 
system with two representatives at each energy 
(Skokos et al. 2002a). As Ej increases, we have one 
branch with z > and the other one symmetric 
with z < 0. Let us consider here the z > case. 
Close to the bifurcating point, the xq initial con- 
dition of the stable xlvl orbit remain close to the 
z = plane. Around the fixed point we find tori, 
which are practically symmetric with respect to the 
z = plane. Away from the bifurcating point, for 
larger values of the energy, the value zq in the ini- 
tial conditions of of xlvl increases. Thus, beyond 
a given energy, in the (x,x,z) projection, the xq 
initial value of the simple unstable periodic orbit 
xl is expected to be away from the tori around 



the stable periodic orbit xlvl, since these tori sur- 
round the xlvl initial condition away from the z=0 
plane. Here we study as an example the case for Ej= 
—5.131377, where the initial conditions of the sim- 
ple unstable periodic orbit xl are (xq, xq, zq, zq) = 
(0.19317510,0,0,0), while these of xlvl are 
lx ,x ,z ,z ) = (0.19221178,0,0.06769306,0). 
Firstly we explore the phase space around the pe- 
riodic orbit xl, which, also at this energy, is simple 
unstable. We do this by increasing the xq initial 
condition by Axo = 0.1, 0.2, 0.3, 0.7 and consid- 
ering the (x, x) surface of section (Fig. [T7|) . The new 
element in Fig. [T7] with respect to the (x,x) phase 
space structure at Ej = —5.1574 (Fig. [7]) is the 
appearance of a chain of 5 islands of stability at 
the location of the third invariant curve. These is- 
lands correspond to the orbit with initial conditions 
(xo + 0.3,0,0,0). Thus, one of the seven invariant 
curves of Fig. [7] has been broken clS 8l result of the 
increase of Ej. We name the rest of the invariant 
curves T la , T 2a , T 4a , T 5a , T 6a ,T 7a , with T la being the 
closest to the initial conditions of xl. 

As in the previous case we investigate the 
neighborhood of the phase space around the 3D 
stable xlvl by applying the same perturbations. 
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Fig. 19. The [x, x) cross-section space at the neighborhood 
of the 3D family xlvl for Ej=-5. 131377. 
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Fig. 20. The (x, x, z) 3D projection of the torus T4 from the 
point view (9,<j>) = (41°, 91°) we consider 10 5 consequents. 

That means we increase the xq initial condition 
of xlvl by Axq = 0.1, 0.2, ... 0.7 and we consider 
the distribution of the consequents in the (x, x, z) 
projection of the 4D space after 10 4 intersections 
with the y=0 plane. We can observe the result in 
Fig. [THJ The red-colored invariant tori around the 
initial condition of the stable p.o. xlvl (red cross 
at (0.19221178,0,0.06769306)) correspond to the 
seven perturbed initial conditions we study. As in 
the case of the perturbed orbits around xl, also for 
xlvl the third set of initial conditions, i.e. the one 
with Axo = 0.3, exhibits a conspicuously different 
dynamical behavior. It forms a set of five small in- 
variant tori, in correspondence to the five islands 
of stability formed on the (x,x) plane in the case 




Fig. 21. 3D projections of torus T5 from the point view 
{0,<t>) = (64°, 127°) (a) and (6,<j)) = (48°, 271°) (b). 

of xl. Each one of the five small invariant tori has 
a rotational torus structure. We name the remain- 
ing six tori that surround the periodic orbit xlvl as 
Ti, T2, T4, T5, Tq, T7. The missing T3 has been sub- 
stituted by the five small invariant tori. In Fig. [18] 
are given also the invariant curves around the p.o. of 
xl, colored green. They lie on the (x,x) plane. We 
can observe that the tori are roughly projected on 
the (x,x) plane in the region occupied by the green 
invariant curves. The projection of the 4D xlvl sur- 
face of section to the (x, x) plane is given in Figfl9l 
We observe that the 2D projection of the tori re- 
sembles the morphology of the invariant curves. An- 
other interesting feature of the T invariant tori at 
the energy Ej= —5.131377 we study now, is that 
the T4 torus to the right of the set of the 5 small 
ones, as well as the tori T5 and Tj are of different 
morphology than the rest of the Tj's. 

In particular the torus T4 (Fig. [18]) has a thin, 
ribbon-like structure. Fig. [20] presents T4 from a 
different point of view {(0,4>) = (41°, 91°)) in the 
(x,x,z) projection of the surface of section and 
helps us understand that T4 indeed intersects it- 
self in this projection at five places. The behavior 
of the T4 is similar with the behavior of the £4 torus 
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at the energy Ej=-5.1574 (Fig. [ID]) . 

T§ has also a thin, complicated, ribbon-like 
structure (Fig. [2Tj) . This time the torus intersects it- 
self only at one region in the (x, x, z) 3D projection 
as we can see in Fig. [2~Ta and in Fig. [2Tb from differ- 
ent points of view of the (x, x, z) projection. Seven 
intersections in the (x,x,z) projection are also ob- 
served in Tj. Ti,T§ and T-j are objects that satisfy 
the definition of the tube tori given by Vrahatis et 
al. (1997). 

As we said, in the (x, x) projection the differ- 
ences in the morphology of T^T^ and Tj from the 
rest Tj's cannot be seen (Fig. [T9|) . Trying to under- 
stand if the difference we observe in the morpholo- 
gies of the 3D projections of the tori reflect some 
morphological differences of the orbits in the config- 
uration space we compare the orbits corresponding 
to the T2 (rotational torus), with those of the tori T4 
and T5 (tube tori) in Fig. [22j Their morphology is 
the expected for quasiperiodic orbits trapped close 
to a stable periodic orbit, in our case xlvl. The 
(x,y) projections of T2,T^ and T5 are practically 
identical with the T^a^T^a and T§ a orbits around xl 
for the same Ej. There is no obvious morphological 
feature that distinguishes the three 3D orbits among 
themselves also in the (x, z) and (y, z) projections. 

T4 offers the opportunity to study in detail 
the structure of one more tube torus in the 4D 
space and compare it with the one depicted in 
Fig. HS1 which we called S4. Applying the color- 
rotation method also in this case, we observe, that 
the structure of T4 (Fig. [23}) is similar with that 
of S4 (Fig. [15J) . We observe in Fig. [23], that mov- 
ing along the tube from A towards B, following the 
directions indicated with arrows, we follow the suc- 
cession of the colors of the color bar from one side to 
the other. At the intersection regions, in the (x, x, z) 
projection, the red color meets blue. This means 
that these regions are not intersections in the 4D 
space. Again here the intersections appear only at 
the 3D projections. 

Different colors at the intersections character- 
ize all tube tori we studied with the color-rotation 
method. A final example is given for the case of T5. 
As we said before the torus T5 is very thin, has al- 
most a ribbon-like structure, and intersects itself at 
one region (Fig. [2Tj) , In Fig. [25] we depict a pro- 
jection in the 3D subspace (x,x,z) and we color it 
according to the values of the consequents in the 4th 
dimension z. By moving counterclockwise from the 
region A we observe again the smooth color varia- 
tion. The colors change from green to light blue and 



then to blue until we reach the region C. At the re- 
gion C T5 intersects itself in the projection (x, x, z). 
Then, moving always counterclockwise, the succes- 
sion of the colors continues as blue — > light blue 
— > green at the region B and finally comes back to 
the region C. At the region C we observe that blue 
color meets red. This means that the points have 
different values in the 4th dimension and C is not a 
region that we have a real intersection of T5 in the 
4D space. If we will continue our counterclockwise 
journey along the thin tubes of T5 we will reach A 
following always a smooth color variation. Of spe- 
cial interest in this case is the folding of the thin, 
ribbon-like torus at the regions A and B. At these 
two regions the internal surface of the torus becomes 
external surface and vice-versa. This warping is en- 
countered along the tubes of several other tori we 
studied. Tj intersects itself in the 3D projections at 
7 regions, while the tubes warp in the journey from 
one intersection to the other. 

4.4.1. Rotation numbers 

The rotation curves in the x-direction along the in- 
variant curves around xl and along the projected 
tori around xlvl in the (x, x) section, as well as 
the rotation curve from the tori around xlvl in 
the (x, x, z) projection, calculated using the defini- 
tion of Eq.(7), are very similar as in the case we 
presented above at Ej = —5.1574. Also here, at 
Ej = —5.131377, there is no special variation of the 
rotation curve at the locations of the tube tori. The 
variation of the rotation curves follows the general 
rules of these curves (see Contopoulos 2002). The 
main difference is that at the location of the sta- 
bility islands around xl and the tori-chains around 
xlvl we find tiny plateaus on the rotation curve as 
expected. 

5. Perturbations in the z- and 
z-direction 

The radial perturbations we studied until now, gave 
us an understanding of the structure of the tori in 
the neighborhood of stable periodic orbits in our 
galactic type Hamiltonian System. However, of par- 
ticular interest is the evolution of the tori also un- 
der vertical perturbations. Especially in 3D galac- 
tic systems we want to know how, and in what ex- 
tent, stars are trapped away from the equatorial 
plane. This study is related with the investigation of 
morphological features as boxy and peanut-shaped 
bulges in disk galaxies (Patsis et al. 2002). 
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Fig. 22. The orbits in the configuration space corresponding to the tori T2 (first row), T4 (second row) and T5 (third row). 
Projections are indicated at the upper left corner of each panel. In the [x, z) and (y, z) projections the scales on the axes are 
not equal, so that we can see the detailed of the orbits. 




5.1. Spaces of section before the 
S — >■ U transition 

Going back to the stability diagram given in Fig. [31 
that describes the stability evolution of the fam- 



ilies xl and xlvl, we study first the phase space 
at an energy before the S —> U transition, i.e. at 
Ej=— 5.207. In this case, on the (x,x) surface of 
section the only 2D simple periodic orbit we have 
with x > is xl and it is stable. Around xl, we find 
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Fig. 24. The (x, x, z, z) 4D surface of section of torus T5 for Ej=— 5.131377. The consequents are colored according to their 
value in the i coordinate. Our view angles are (6, <j>) = (30°, 72°). 



invariant curves as expected, if we perturb its xq ini- 
tial condition. As an example we consider the invari- 
ant curve, which results from the initial conditions 
(x + Ax ,x ,z ,zo) = (0.18312784 + Ax , 0,0,0), 
with Axo=0.1. If we "perturb" this quasi-periodic 
orbit, i.e. if we increase further the zq = initial 
condition by Azo=0.1, 0.2 . . . , and integrate again 
our initial conditions, we encounter tori around xl. 
These are both rotational and tube tori. The tori, 
for small Azo surround the invariant curve we per- 
turbed. For Azo >0.3 the tori appear distorted and 
inclined with respect to the z = plane. However, 
for Azo=0.41 the dynamical behavior at the neigh- 
borhood of the quasi-periodic orbit changes. About 
the first 900 consequents of this orbit form an ob- 
ject that resembles a rotational torus. By this we 
mean that initially the consequents stay approxi- 
mately on a toroidal surface. However, the following 
consequents diffuse in the phase space and occupy 
a larger volume in it. In Fig. [25a we observe the 
(x,x,z) projection of the 4D surface of section for 
this orbit, where the 4th dimension is represented 
by the colors of the 4th i coordinate. The num- 
ber of consequents in Fig. [25b , is 1100. They form 
a (rotational) toroidal surface, with a smooth color 
variation around it as we can see by comparing the 
colors on its surface with the color bar to the right 



of the figure. We note, that 900 consequents do not 
suffice to fill densely the surface of the toroidal ob- 
ject. Thus, on Fig. [25k we observe distinct points 
rather than a toroidal object. However, starting e.g. 
from upper left, we observe the succession of colors 
from green — > yellow — > orange — > red — > orange 
— > yellow — > green — > to blue shades. In Fig. l25a 
some points have already started deviating from the 
"torus" and soon they occupy a larger area of the 
phase space. In Fig. [25b . the dense cluster of red 
points corresponds to the 900 consequents we plot- 
ted colored in Fig. l25h . The rest of the red points 
that occupy a larger area of the phase space be- 
long to the same orbit, which has been integrated 
now for 1.2 x 10 4 consequents. In the same figure 
we give as well the invariant curve around xl at 
this energy, drawn with green color. This dynamical 
behavior indicates stickiness (Contopoulos & Har- 
soula 2008). Qualitatively, we find similar results if 
we increase zq instead of zq, in the initial condi- 
tions (x ,x ,z ,zo) = (0.28312784,0,0,0) following 
the procedure described above. 
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Fig. 25. The 3D (x,z,x) projection of the 4D surface of section in the neighborhood of xl for Ej =— 5.207. The initial 
conditions of the orbit are (xq, xq, Zq, Zq) = (0.28312784, 0, 0.41, 0) (see text), (a) the first 1100 consequents colored according 
to their position in the 4th dimension z. The point of view is (9, (f>) = (78°, 82°). (b) The first 1.2 xlO 4 consequents (red dots). 
The initial condition of xl is indicated with a green dot, while the green invariant curve around it correspond to the orbit with 
initial conditions (0.28312784,0,0,0). The red points diffuse in the phase space after staying on a toroidal surface for about 
900 consequents. 



5.2. Spaces of section after the 
S — >■ U transition 

As we already have seen, after the S — > U tran- 
sition, xl becomes simple unstable, while the 3D 
stable family xlvl has been bifurcated. The en- 
ergy Ej=— 5.131377 is already beyond the bifurcat- 
ing point (Fig. [3j). We have examined the phase 
space structure in the neighborhood of xl and xlvl 
by perturbing their xq initial conditions and the 
result is depicted in Fig. [181 At energies not far 
away from the bifurcating point A (Fig. [3]), as is 
Ej=— 5.131377, we could always find a "±Azo" per- 
turbation of the quasi-periodic orbits, so that the 
initial condition (xq + Axo, 0, ±Azq,0) gives a torus 
surrounding the invariant curve (xq + Axq, 0,0,0), 
with (xq, 0,0,0) being the initial condition of xl. 



However, at larger energies this was not always 
possible. We could find invariant curves around xl, 
that deviate in the phase space from the invariant 
tori around xlvl at the same energy. For example 
at Ej=— 4.98996, we find invariant curves around 
xl and invariant tori around xlvl by perturbing 
the xq initial conditions of both families. In Fig. [26] 
we give in the 3D (x,z,x) projection the invariant 
curves around xl for xq + Axq = 0.1, 0.2 and 0.3 
(green curves) and the corresponding tori around 
xlvl applying the same perturbation to its initial 
conditions. The invariant tori "float" above the in- 
variant curves. We observe that as we depart from 
the xlvl initial condition the invariant tori become 
thicker in z. Beyond a certain xq + Axq perturba- 
tion there are invariant tori that intersect the (x, x) 
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Fig. 26. The 3D (x,z,x) projection of the 4D surface of section at Ej=— 4.98996 for three invariant tori around xlvl and 
three invariant curves around xl (green lines). Our point of view is (9, <j>) = (56°, 280°). The initial condition of xl is indicated 
with a green star and the initial condition of xlvl is given with a red "x". 




I.b D.S ,.S C.5 I !.S 

\ 1 

Fig. 27. (a) A chaotic orbit at Ej=-4.98996 in the 3D (as, z, x) projection (point of view (9, tf>) = (55°, 263°)) of the 4D surface 
of section is given with blue dots. We also give three rotational tori around xlvl (red), and seven invariant curves around xl 
in the (as, x) plane, (b) The 4D surface of section of the cloud of blue points we give in (a) (point of view (9, <f)) — (60°, 45°)). 
It is characterized by scattering of their colors in the 4th dimension, (c) The (x,x) surface of section for the orbits in (a). We 
see that the blue points extend inside as well as outside the innermost rotational torus, (d) The chaotic orbit (blue points) 
together with the quasi-periodic xlvl orbit corresponding to the innermost rotational torus in (a). 



plane (not depicted in Fig. [26]) . 

In the case depicted in Fig. [26] we succes- 
sively "perturb" the initial conditions of the quasi- 



periodic orbits that correspond to the red invariant 
tori by — Azq trying to find other invariant tori, that 
their projections in z reach the (x, x) plane. Follow- 
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Fig. 28. (a) The (y, z) profile of the orbit corresponding to the cloud of blue points in Fig. 127b . (b) The same orbit in the 
configuration space from a point of view (8,<j>) — (63°, 280°). 



ing this procedure we find a last invariant torus with 
(s + 0.1, 0, z + Az , 0) = (0.3160689, 0, -0.08, 0), 
i.e. for Azq = —0.08. For Azq = —0.09 instead of a 
torus we find a cloud of points at the area between 
the invariant curves and the invariant tori. This 
cloud is given in Fig. [27b with 1.2x10 blue points, 
together with the three rotational tori around xlvl 
depicted in Fig. [26] (colored red) and seven invariant 
curves around xl on the (x,x) plane (black curves). 
For about the first 3000 consequents, the cloud has 
a vague toroidal structure in the (x, x, z) projection. 
In Fig. 127b we apply the color-rotation method in 
the set of blue points of Fig. [27k . which we color 
according to their value in the z coordinate. It be- 
comes clear that the points do not have a smooth 
variation in the 4th dimension. The same color mix- 
ing is also present in the figure with fewer conse- 
quents. Particularly helpful is the (x, x) projection 
of Fig. l2Th . given in Fig. [2Tb . In this projection we 



can clearly observe, that blue points of the cloud 
can be found inside, as well as outside of the red 
innermost invariant rotational torus. The vast ma- 
jority of the points diffuses outside the the (x, x) 
projection of that torus between 4xl0 3 to 5xl0 3 
consequents, building a second ring at the outer 
side of the torus of lower intensity than the ring in 
its inner side. It is interesting to realize that such 
chaotic orbits may contribute to the thickening of 
observed structures as the peanut-shaped bulges of 
disk galaxies. The backbone of such peanut-shaped 
structures is the xlvl family (Patsis et al. 2002). 
In Figj27U the trapped quasi-periodic orbit around 
the stable p.o. xlvl at this energy (Ej=— 4.98996) 
is plotted with red dots, while the blue points give 
the chaotic orbit that we plot in Fig. 127b in the con- 
figuration space ((x,z) projection). The blue orbit 
in FigJ27l is integrated for a time corresponding to 
0.2 Hubble times, i.e. in a time generally believed 
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Fig. 29. The relative location of the consequents of the chaotic orbit of Fig. 1281 with respect to the location of the tube torus 
"X4" (see text). In (a) and (b) we give with a blue cross, the location of a single consequent that it is projected in the interior 
of the (x,x) and (z, z) projections of "X4" respectively. In (c) and (d) we give all the 1.2 x 10 4 consequents of the orbit, clearly 
showing that we have points inside and outside "X4" . 



that morphological features as the galactic bars can 
well survive (see e.g. Debattista et al. 2006). 

For longer times, of the order of 0.7 of a Hub- 
ble Time, the shape of the "blue" orbit still fol- 
lows the peanut morphology, having populated both 
branches of the stable p.o. xlvl (Patsis et al. 2002) 
at this energy. Thus, the peanut shape is fully devel- 
oped, despite the fact that the orbit reaches z values 
1.5 times larger than the quasi periodic orbit in the 
neighborhood of xlvl (innermost rotational torus). 
In Fig. [28b we give the (y, z) projection for compari- 
son with the models of periodic orbits in Patsis et al. 
(2002), while in Fig. 128b we give the same orbit from 
the point of view ((0,0) = (63°, 280°)), so that it 
becomes clear the relation between the peanut side- 
on profile of the orbit and its "pencil-sharpener" 3D 
morphology. This class of orbits may be very impor- 
tant for the kinematics of peanut-shaped bulges of 
disk galaxies, as they increase the dispersion of the 
velocities of the stars, participating in this struc- 
ture without destroying its morphological profile. 
These orbits are chaotic, but for times important 
for Galactic Dynamics have consequents that stay 



close to the invariant tori of xlvl in the phase space. 
However, for larger times they diffuse away from the 
invariant tori in the phase space. 

None of the blue consequents presented in the 
(x,x) projection (Fig. 127b) inside the innermost red 
torus could be found projected simultaneously in 
the interior of the torus in the (z, z) projection. 
However, this was not true for the location of con- 
sequents relative to tori at larger distances from the 
periodic orbits than the three rotational tori we con- 
sidered in Figs. [26] and [271 An example is given for 
the torus we obtained by perturbing by Axo=0.4 
the xlvl initial conditions. This is a tube torus and 
let us call it "XT' . In Fig. [29] we give with a blue 
cross the position of a consequent which is simul- 
taneously projected in the interior of "XT' of the 
(x, x) (Fig. [29k ) as well as of the (z, z) projection 
(Fig. [29b). It is one of the consequents depicted in 
Figs. [29b and Fig. [29H . where we clearly observe 
that we have consequents of the "blue" orbit in- 
side and outside the tube torus "XT'. This tells 
us that in this case takes place the phenomenon 
of Arnold diffusion (see e.g. Contopoulos 2002, pg. 
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Fig. 30. Stability Diagram for -4.5 <Ej < -4.37, that gives 
the stability indices of the family xl and its bifurcating family 
of p.o. xlvl, in the axisymmetric case of our model (q a = 
1,96 = 1)- 



344). Nevertheless the time scale for this exceeds 
one Hubble Time. Thus, in this particular case the 
consequences for the dynamics in the neighborhood 
of xlvl are negligible. However, in a forthcoming 
paper we investigate in detail this as well as other 
cases of Arnold diffusion in the Hamiltonian system 
of Eq. 4 and Eq. 5. 

Before closing this section we want to underline 
that the same procedure can be repeated by per- 
turbing the zq initial condition, reaching the same 
results as regards the appearance and the foliation 
of the tori in the 4D space. 

6. The evolution of the tube tori as 
the perturbation increases 

Finally we examine the evolution of the invariant 
tori as the perturbation of our system increases, 
starting from the axisymmetric case (q a = 1 and 
qb = 1). In our model we increase the perturbation 
by increasing the triaxiality of the system through 
the parameters q a and qb- At first we examine the 
axisymmetric case and then we introduce a small 
perturbation (q a = 1.01 and qb = 1), in order to 
study possible qualitative differences that are in- 
troduced in the system by the perturbation. 

6.1. The Axisymmetric case 
(q a = 1, q h = 1) 

In 3D axisymmetric potentials the stability index 
of the central family, which is related with pertur- 




Fig. 31. (x, x, z) 3D projections of the 4D surface of section 
in the axisymmetric case, for Ej=— 4.4461, at the neighbor- 
hood of the marginally stable p.o. xlvl. (a) The first 8 invari- 
ant tori VI. .. V8 (point of view (0,4>) = (38°, 211°)). (b) The 
invariant tori V9. . . V15, (point of view (9, (j>) = (36°, 244°)). 
(c) The last six tori of our sample, V16. . . V21, (point of view 
(0, <t>) = (49°, 195°)). We observe the deformation of the last 
torus V21. 



bations vertical to the equatorial plane, becomes 
tangent to the b=— 2 axis at the vertical resonances 
(e.g. Patsis & Grosb0l 1996). Fig. [301 gives the evolu- 
tion of the stability of the central family of periodic 
orbits of our system, xl, and its vertical bifurcation 
(xlvl) for -4.5 <Ej< -4.37. We observe that xl 
(bl, b2 indices with black lines) is initially stable 
and at Ej=— 4.4472 we have a tangency of bl with 
the —2 axis. This tangency corresponds to the ver- 
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Fig. 32. The (x, x, z) 3D projection of the 4D surface of 
section in the neighborhood of the p.o. xlvl for Ej=— 4.39. 
Our point of view is (0, (p) = (27°, 316°). 



tical 2/1 resonance. At this point it is bifurcated a 
new family (xlvl) (the two stability indices of xlvl 
are drawn with red lines), which remains marginally 
stable, since we find that one of its stability indices 
stays on the —2 axis. 

We have taken several surfaces of section and 
we studied the structure of the phase space at the 
neighborhood of the marginally stable p.o. xlvl. 
We first studied the phase space structure close to 
xlvl for Ej=— 4.4461, just after the tangency of 
the bl index with the b=— 2 axis. We integrated 
21 initial conditions at the neighborhood of xlvl 
(x ,z ,xo,z ) fa (1.298965,0.041883918, 
-.0036186582,0.026639957), by perturbing the x 
initial condition of xlvl by Ax = 0.02, 0.04 . . . 0.42 
successively. Fig. [31] gives the (x, x, z) 3D pro- 
jections of the tori we found around xlvl. In 
Fig. [31"a we observe the first eight rotational tori 
Vi,V2,...Vs and in Fig. [31b the rotational tori 
Vg . . . V15 . Finally in [3Tb we give the last six tori 
Vi6 ■ ■ ■ V20 and V21, of our sample. Each torus con- 
sists of 10 4 consequents. In this dense coverage of 
the Axo = 0.4 space away from the xq initial con- 
dition of the xlvl p.o., we encountered only rota- 
tional tori nested around xlvl, until we reached 
V21. There, perturbing xq by Axo = 0.42, the mor- 
phology of the torus changes. The V21 torus has a 
larger projection on the z-axis and does not follow 
anymore the pattern of the nested tori. 



We repeat the same procedure at a larger en- 
ergy, Ej=— 4.39. This time the x initial condition of 
xlvl is perturbed successively by Ax = 0.1, 0.2, 0.3 
and 0.4. The result is given in Fig. [32l having a 
point of view (9,<f>) = (27°, 316°). Four rotational 
tori surround the xlvl p.o. and are indicated with 
arrows (Ui,U2,Us and U4). Each torus consists of 
10 4 consequents. The tori are distorted, with the 
last one having a morphology similar to V21 of the 
previous case. 
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Fig. 33. Stability Diagram for -4.7 < Ej < -4.3, that 
shows the stability of the family xl and its bifurcating family 
of p.o. xlvl. 



6.2. The small perturbation case 
q a = 1.01,g 6 = 1.0 

In order to study qualitative differences from the 
axisymmetric case when a small perturbation is in- 
troduced in the system, we increase the q a param- 
eter from q a = 1 to q a = 1.01. Fig. [33] gives for this 
system the evolution of the stability of the central 
family xl and that of its bifurcations for —4.7 < Ej 
< —4.3. We observe that xl (bl, b2 indices with 
black lines) is initially stable and at Ej=— 4.54 we 
have a transition from stability to simple instability. 
The family xl becomes again stable at Ej=— 4.386. 
At the transition from stability to simple instabil- 
ity xlvl is bifurcated (red lines show its stability 
indices ), as stable. At a larger energy (Ej=— 4.525) 
we have a S — > A transition to complex instability 
of the family xlvl. In this case the stability indices 
become complex numbers and they do not appear 
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W5 



Fig. 34. The (x,x,z) 3D projection of the 4D surface of section for Ej=— 4.4461 close to the stable p.o. xlvl. Our point of 
view is (6,(f>) = (55°, 194°). 



in the figure. There is a narrow complex unstable 
interval of xlvl until Ej=— 4.525, where xlvl be- 
comes again stable. In the present paper we study 
the phase space only in the neighborhood of stable 
p.o. 

We choose a value of energy for which we have 
studied the orbital dynamics at the neighborhood 
of a stable xlvl orbit at the axisymmetric case, 
Ej=— 4.4461, and we perturb first he xq initial con- 
dition of xlvl by Ax = 0.02,0.04 ...0.14 succes- 
sively. 

In Fie. 1341 we observe seven tori at the ( X • X - Z J 
projection, which we name W\,W2,---Wj, sur- 
rounding the fixed point of the xlvl p.o. at this Ej. 
Each torus consists of 10 4 consequents. The mor- 
phology of Wi,W2,Wz,Wi and Wq indicates that 
they are rotational while the morphology of W§ and 
W-j is typical for tube tori. The tubes of W5 have 
many self-intersections. However, the gaps between 
them do not fill even if we consider 10 6 consequents. 
Thus, we conclude that W5 is a tube torus. On 
the other hand, despite its wavy character, W7 is 
a tube torus with only one self intersection in the 
3D (x,x,z) projection. We note that tube tori ap- 
pear in the 3D projections of the spaces of section 
as soon as a perturbation is introduced, even if it is 
a small one. 

7. Summary and Conclusions 

In this phenomenological paper we have studied in 
detail the structure of phase space in the neighbor- 
hood of stable periodic orbits in a 3D potential that 



represents a rotating, thick galactic disk. We have 
visualized the 4D spaces of section by means of the 
color and rotation method. This allowed us to clar- 
ify the properties of the invariant tori that we en- 
counter in the vicinity of the stable p.o. We have 
also examined chaotic zones that we have found 
encaged between tori and the effect that these zones 
could have in the support of structures observed in 
thick galactic disks. Below we summarize our con- 
clusions and compare our results with the results of 
previous studies. 

(1) On the equatorial plane (z=0), the dynamics of 
the system is determined by the presence of the 
central family xl. Even in the intervals that xl 
is characterized as simple unstable (U), we find 
in the (x,x) surface of section invariant curves 
surrounding the periodic orbit. 

(2) In general we find around stable periodic or- 
bits invariant tori. In the 3D projections we 
encounter the invariant tori that have been 
found in the papers by Froeschle (1970, 1972), 
Martinet and Magnenat (1981), Contopoulos et 
al. (1982), Magnenat (1982), Patsis & Zachi- 
las (1994) and Vrahatis et al. (1997), named by 
the latter authors "rotational tori". However, 
we have found as well the other kind of invari- 
ant tori, named by Vrahatis et al. (1997) "tube 
tori" . The fact that this second kind of tori ap- 
pear frequently in the vicinity of stable periodic 
orbits in totally different physical systems as are 
Hamiltonian systems of galactic type (this pa- 
per) and 4D symplectic maps associated with 
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the problem of beam stability in circular parti- 
cle accelerators (Vrahatis et al. 1997), indicates 
that they represent a generic behavior in 3D 
orbits. 

(3) The method of color and rotation was particu- 
larly enlightening in understanding the proper- 
ties of the two kinds of tori. 
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• As regards the rotational tori in all cases we 
studied in our potential the continuation of 
color, that guarantees the smooth variation of 
the 4th coordinate, changes side along certain 

lines (Figs. QUQ2D. This shows that in the 3D 
projections we have intersections in the rota- 
tional tori, which, in all cases we encountered 
in this study, appear four times on the tori. 

• As regards the tube tori, their characteristic 
feature, the self-intersections of their tubes, was 
a projection effect in all cases, since at the in- 
tersections meet always different colors. We re- 
mark that we found always an odd number of 
self-intersections in the 3D projections of the 
tube tori. In the J^D space we did not find any 
self-intersections of the tube tori. 



We found that the two kinds of tori cannot be 
distinguished neither from their rotation num- 
bers, nor from the morphological features of the 
orbits in the configuration space. 

(4) As we depart from the initial conditions of the 
stable p.o., the 3D projections of the tori ap- 
pear distorted. Beyond a certain torus we found 
sticky chaotic orbits that for a long time remain 
on toroidal surfaces on which the consequents 
have a smooth color variation and finally dif- 
fuse in phase space, where we observe mixing 
of colors. 

(5) We found chaotic orbits, which can be encaged 
by tori in the 4D space for times longer than a 
Hubble time, but not forever, since there are no 
topological barriers. We presented an example, 
where such orbits support a peanut morphol- 
ogy, similar but thicker than the morphology 
supported by the corresponding quasi-periodic 
orbits trapped around the stable xlvl. 

(6) Tube-tori have been found only in the 3D pro- 
jections of orbits in perturbed systems. In the 
axisymmetric case we encountered only rota- 
tional tori. 
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